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Abstract 

A new 3.2-3.5 /im spectrum of the protostar Mon R2/IRS-3 confirms our previous 
tentative detection of a new absorption feature near 3.25 /im. The feature in our new 
spectrum has a central wavelength of 3.256 fxm (3071 cur 1 ) and has a full-width at half 
maximum of 0.079 /im (75 cm -1 ). We explore a possible identification with aromatic 
hydrocarbons at low temperatures, which absorb at a similar wavelength. If the feature 
is due to aromatics, the derived column density of C-H bonds is ~1.8 x 10 18 cm -2 . If 
the absorbing aromatic molecules are of roughly the same size as those responsible for 
aromatic emission features in the interstellar medium, then we estimate that ~9% of the 
cosmic abundance of carbon along this line of sight would be in aromatic hydrocarbons, in 
agreement with abundance estimates from emission features. 



Subject headings: infrared: general — ISM: dust, extinction — ISM: molecules 
— stars: pre-main sequence 
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1. Introduction 



The C-H stretch absorptions of many of the organic molecules expected to be formed or 
condensed on molecular cloud dust lie in the 3.2-3.6 /im region, on the long wavelength side 
of the 3.1 fim H 2 ice band which dominates the spectrum of embedded sources. Sellgren, 
Smith, & Brooke (1994) recently reported a tentative detection of a new absorption feature 
at 3.25 /im (3078 cm^ 1 ) toward Mon R2/IRS-3, a protostar in the Mon R2 star formation 
region (Beckwith et al. 1976). Their spectrum had a resolution A/AA ~ 720 at 3.25 /jm. 
Here, we present a new spectra of Mon R2/IRS-3 with a resolution of 1000 which confirms 
the presence of a 3.25 /im feature. Some possible identifications are discussed. 



2. Observations 



The latest observations of Mon R2/IRS-3 were made on 1994 October 8 at the United 
Kingdom Infrared Telescope (UKIRT) on Mauna Kea. The CGS4 long-slit spectrometer 
(Mountain et al. 1990) was used with the 75 lines mm -1 grating in first order and the 
300-mm focal length camera. This provided a wavelength resolution of 0.0033 /im (A/AA 
= 1000 at 3.25 /im). The spectrometer is designed to have only one resolution element 
per pixel, so improved sampling of the spectrum was obtained by moving the detector 
by one-third of a resolution element between individual spectra and repeating this until 
two resolution elements were observed by each pixel. The observations consist of two 
overlapping grating positions, at 3.16-3.37 /tm and 3.34-3.55 /xm. The pixel size was 1.55". 
The spectrometer slit was 90" x 1.55" with the long direction oriented east-west. The 
sources were nodded ~12" along the slit for background subtraction. An argon spectrum 
in second order was used for wavelength calibration. We compared our spectrum of Mon 
R2/IRS-3 with the star HR 1948 (09Iab:) for atmospheric cancellation. The airmass 
difference between Mon R2/IRS-3 and HR 1948 was always less than 0.03. 

In the final spectra, several points at 3.313 - 3.321 /im affected by strong telluric CH 4 
have been removed. We have also removed points near 3.297 fim which may have been 
affected by any photospheric Pfund S feature in the 09Iab: atmospheric comparison star. 
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3. Results 



The new spectrum of Mon R2/IRS-3 is shown in Figure 1. The observations fall 
in the region of the 3.1 /im H2O ice band and the broad absorption wing which peaks 
near 3.3-3.4 /im (Smith, Sellgren, & Tokunaga 1989). The intrinsic spectral shape of this 
absorption is uncertain. Thus the best continuum to use for deriving the optical depth 
of narrow absorption features in this region is a local continuum which passes smoothly 
through those parts of the spectrum not containing narrow absorption features. We have fit 
a second-order polynomial to the spectrum of Mon R2/IRS-3, excluding data at 3.2-3.3 /im 
and longward of 3.4 /im from the fit. The choice of excluded regions is the same as that 
used by Sellgren et al. (1994). Our adopted continuum is shown as a solid line in Figure 1. 

The derived optical depth is also shown in Figure 1. We fit two Gaussians to the optical 
depth curve. The central wavelength, full width at half-maximum (FWHM), and optical 
depth of each Gaussian were varied to produce the best fit to our observations. We derive 
central wavelengths of 3.256 ± 0.003 fim and 3.484 ± 0.003 /im (3071 ± 3 cm" 1 and 2870 
± 2 cm -1 ) for the 3.25 /im and 3.48 /im features, respectively. We also find FWHM values 
of 0.079 ± 0.007 and 0.117 ± 0.007 /im (75 ± 6 cm" 1 and 97 ± 6 cm' 1 ) for the 3.25 
/im and 3.48 /im features, respectively. Our new measurements of the central wavelengths 
and widths agree well with those of Sellgren et al. (1994). The 3.25 /im optical depth we 
measure, 0.045, also agrees well with Sellgren et al. (1994). The 3.48 /im optical depth we 
derive, 0.058, does not agree with the value of 0.036 measured by Sellgren et al. (1994). 
However, the optical depth is sensitive to the choice of continuum, so the Sellgren et al. 
(1994) spectrum provides the most reliable value for the 3.48 /im optical depth because 
the current spectrum (Fig. 1) does not extend to long enough wavelengths to provide 
continuum on the long wavelength side of the 3.48 /im feature. 



4. Discussion 



The 3.48 /im feature was first identified by Allamandola et al. (1992) toward four 
protostars. They attributed the feature to C-H bonds in hydrocarbons with "diamond-like' 
bonding. This feature in Mon R2/IRS-3 and other sources is discussed in more detail by 
Brooke, Sellgren, & Smith (1995). 
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Standard references on room temperature infrared spectra suggest that the 3.25 /zm 
feature might be due to a C-H stretch of the =CH 2 group in an alkene, which occurs at 
3.23-3.25 /iin (e.g. Williams & Fleming 1987). An alkene identification, however, is unlikely 
because alkenes have a second, comparably strong, feature at 3.29-3.32 /zm which is not 
observed toward Mon R2/IRS-3. We have searched the low temperature laboratory spectra 
of pure ices and ice mixtures with compositions thought to be appropriate to molecular 
clouds (d'Hendecourt & Allamandola 1986; Grim et al. 1989; Hudgins et al. 1993). These 
spectra reveal no obvious absorption features near 3.25 /zm. 

We suggested earlier (Sellgren et al. 1994) that the 3.25 /zm feature may be due 
to absorption by aromatic hydrocarbons at low temperature, based on a similarity in 
wavelength to the C-H stretch of polycyclic aromatic hydrocarbons (PAHs) isolated in 
neon matrices at a temperature of 4.2 K (Joblin et al. 1994). The aromatic C-H stretch 
wavelength is a function of temperature, increasing with increasing temperature (Colangeli, 
Mennella, & Bussoletti 1992; Joblin et al. 1994, 1995). Aromatic hydrocarbons are a 
promising candidate for the 3.25 /im absorption feature, since aromatic emission features 
at 3.3, 6.2, 7.7, 8.6, and 11.3 fim have been observed throughout the interstellar medium 
of our own and other galaxies. Corresponding absorption features have been searched for, 
but until now have not been definitely detected in molecular clouds. The infrared emission 
features have been attributed to a variety of aromatic substances, including hydrogenated 
amorphous carbon (HAC) grains (Blanco, Bussoletti, & Colangeli 1988; Ogmen & Duley 
1988), PAHs (Leger & Puget 1984; Allamandola, Tielens, & Barker 1985), quenched 
carbonaceous composite (QCC) grains (Sakata et al. 1987), and other aromatic materials 
(see Sellgren 1994 for a review of proposed identifications). 

We compare in Table 1 the observed wavelength of the 3.25 fim feature toward Mon 
R2/IRS-3 with the wavelengths of several aromatic substances. We list in Table 1 the 
measured wavelengths of solid QCC (Sakata et al. 1990), the 3.3 /im interstellar aromatic 
emission feature (Tokunaga et al. 1991), the PAH molecule coronene in the condensed 
phase and the gas-phase (Flickinger, Wdowiak, & Gomez 1991), solid HAC (Biener et al. 
1994), and the PAH molecules coronene and pyrene isolated in a neon matrix (Joblin et al. 
1994). 

Joblin et al. (1995) have examined the temperature dependence of the C-H stretch 
wavelength of gas-phase aromatic molecules in detail. They state that the wavelength 
increases with increasing temperature due to anharmonic coupling of the C-H stretch mode 
with excited longer wavelength modes. In Table 1 we also present the predicted wavelengths 
for each aromatic material, when shifted from the temperature at which the measurement 
was made to a temperature of 80 K, appropriate for the icy grains toward Mon R2/IRS-3 
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(Smith et al. 1989), using Eq. 5 of Joblin et al. (1995) and the assumption that the 
neon matrix does not introduce a wavelength shift from the gas phase. The temperature 
dependence of the aromatic C-H stretch wavelength (Joblin et al. 1995) was derived for 
gas-phase aromatic molecules, and we caution that solid-phase aromatics, such as HAC or 
QCC, may not follow the same relation. 

In Figure 1, we compare the optical depth profile of the 3.25 /iin absorption feature 
and the profile of the 3.3 /mi aromatic interstellar emission feature in IRAS 21282+5050 
(Nagata et al. 1988), after continuum subtraction (Tokunaga et al. 1991), and after shifting 
the center of the emission feature to the predicted wavelength at 80 K (see Table 1). The 
two feature profiles show reasonable agreement, although since the width of each feature is 
probably dominated by different processes, such agreement may be fortuitous. 

The average of the observed feature wavelengths from this paper and Sellgren et al. 
(1994) is 3.253 ± 0.004 /mi, which is shorter than the aromatic hydrocarbon wavelengths 
in Table 1 by 0.004-0.032 /mi. The fact that the 3.25 /tin absorption feature just barely 
overlaps the short wavelength side of the range of cold aromatic hydrocarbon wavelengths 
presents a problem, since moving the aromatic C-H vibration to shorter wavelengths (higher 
frequencies) means strengthening the C-H bond, something that seems difficult to achieve 
if the aromatic hydrocarbons are immersed in an ice matrix of some sort. 

Any identification of the 3.25 /mi feature at this time rests only on one absorption 
feature, and the wavelength match with aromatic hydrocarbons is not exact. A search for 
the longer wavelength features associated with aromatic hydrocarbons would provide one 
test of this identification. 

If we assume that the 3.25 /iin absorption feature is due to aromatic hydrocarbons, 
the column density of aromatic C-H bonds along the line of sight to Mon R2/IRS-3 can be 
estimated. Measurements of aromatic hydrocarbons in absorption are important because 
estimates of the abundance of aromatic hydrocarbons from the observed emission features 
(Allamandola et al. 1989; Puget & Leger 1989; Joblin, Leger, & Martin 1992) are much less 
straightforward. 

To estimate the column density of aromatic C-H bonds, we use the relation, N ~ 
tAu/A, where r is the maximum optical depth of the 3.25 /tm absorption feature, Av is 
the feature FWHM in cm -1 , A is the integrated absorbance, and N is the derived column 
density of molecular bonds (Allamandola et al. 1992). An average of the results of this 
paper and Sellgren et al. (1994) gives r(3.25 /im) = 0.047 and Av = 66 cm -1 for the 3.25 
/im feature. For the three aromatic molecules, pyrene, coronene, and ovalene, studied by 
Joblin et al. (1994), the value of A per aromatic C-H bond for the 3.25 /tm feature was 
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0.7-1.4 x 10~ 18 cm bond -1 in the solid phase and 2.1-4.1 x 10~ 18 cm bond -1 in the gas 
phase. We average over all three molecules in both phases, to estimate an average value of 
A = 1.7 x 10~ 18 cm bond -1 . We thus derive a column density of aromatic C-H bonds of 
iV(C-H) ~ 1.8 x 10 18 bonds cm -2 along the line-of-sight. 

The abundance by number of aromatic C-H bonds, X(C-H), is the ratio of the column 
density of aromatic C-H bonds divided by the total hydrogen column density, N H . We 
estimate Nh in two ways. The silicate optical depth, r(9.7 /jm) = 4.3, observed toward 
Mon R2/IRS-3 (Willner et al. 1982) implies Ay = 80 mag assuming A y /r(9.7 //m) = 18.5 
(Mathis 1990). However, Ay/r(9.7 /im) is observed to vary by a factor of two (Mathis 
1990). An independent estimate of Ay comes from the 4.6 /im 13 CO gas absorption 
observed toward Mon R2/IRS-3 (Mitchell 1995), which gives iV( 13 CO) = 1.6 x 10 17 cm -2 . 
If we assume A v /N( 13 CO) = 4 x 10~ 16 cm 2 mag (Dickman 1978), then the 13 CO gas 
column density implies Ay = 64 for Mon R2/IRS-3, in good agreement with the value 
derived from the silicate feature. We then convert our average value of Ay = 72 to Nh 
by assuming N H /A V = 1.9 x 10 21 cm -2 mag -1 (Mathis 1990). This implies N H = 1.4 x 
10 23 cm -2 for Mon R2/IRS-3. Again there is some uncertainty in this because the value of 
Nh/ Ay measured in the diffuse interstellar medium may not hold in molecular clouds. Our 
derived value of Nh implies that X(C-H) = 1.3 x 10~ 5 toward Mon R2/IRS-3. For a solar 
abundance of carbon, X(C)/X(H) = 3.6 x 10~ 4 by number (Anders & Grevesse 1989), our 
estimate of X(C-H) implies that ~ 4% of the total carbon along the line of sight toward 
Mon R2/IRS-3 is locked in aromatic C-H bonds. 

The total number of carbon atoms in aromatic hydrocarbons will be larger. If the 
absorbing aromatic hydrocarbons have the same size distribution as the emitting aromatic 
hydrocarbons, then we can use model results for the interstellar emission features to 
estimate the fraction, /, of the number of carbon atoms in aromatic C-H bonds, compared 
to the total number of aromatic carbon atoms. The value of / depends on the aromatic 
hydrocarbon size, with a smaller value for larger aromatic hydrocarbons. Desert, Boulanger, 
& Puget (1990) present a model of interstellar dust, including size distributions for different 
grain components and an analytic approximation for / as a function of radius a for 
PAH molecules. We have used their model, with a = 4-12 A for PAHs, to calculate a 
size-averaged value for / of 0.40. The value of / for the absorbing aromatic hydrocarbons 
also depends on the degree of dehydrogenation in the interstellar medium, but aromatic 
hydrocarbons are predicted to be fully hydrogenated in molecular clouds shielded from 
ultraviolet radiation (Allamandola, Tielens, & Barker 1989). Thus the total amount of 
carbon in aromatic hydrocarbons is roughly a factor of ~2.5 times higher than the amount 
of carbon participating in aromatic C-H bonds. If the 3.25 /im feature is due to absorbing 
aromatic hydrocarbons with a size distribution similar to that adopted by Desert et al. 
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(1990) for the emitting aromatic hydrocarbons in the interstellar medium, this would make 
the fraction of carbon in aromatic hydrocarbons ~9%. If the absorbing aromatic C-H bonds 
are instead attached to larger structures, for instance if the aromatic absorption is due to 
hydrogen on the surfaces of large amorphous carbon grains while the aromatic emission is 
due to small PAH molecules, then the fraction of carbon in such structures would be much 
larger than we estimate from the Desert et al. (1990) model. 

Our estimate of the carbon abundance in aromatic hydrocarbons of ~9% falls within 
the range of previous estimates for the aromatic hydrocarbon abundance, which vary from 
0.8% to 18% of the total carbon abundance (Lepp et al. 1988; Allamandola et al. 1989; 
Puget & Leger 1989; Joblin, Leger, & Martin 1992). Thus if the 3.25 /im feature is due 
to aromatic hydrocarbons, we estimate that a significant fraction of carbon remains in 
aromatic hydrocarbons in molecular cloud dust. 

If the 3.25 /im feature is due to, or contains contributions from, non-aromatic species, 
then the abundances of aromatic hydrocarbons along the line-of-sight derived above become 
upper limits. If it can be shown that none of the feature is due to aromatic hydrocarbons, 
then the abundance of carbon trapped in aromatic hydrocarbon molecules may be much 
lower in molecular clouds than in photodissociation regions or the diffuse interstellar 
medium. Aggregation of aromatic hydrocarbon molecules into larger graphitic-like 
structures is one possible explanation. 

The most pressing need is to detect the 3.25 /im feature in other sources, both 
protostars and field stars behind molecular clouds. Brooke et al. (1995) have recently 
detected the 3.25 /im feature toward the protostars NGC 7538/IRS-l and S 140/IRS-l, 
but observations are needed over a wider range of physical conditions. This will determine 
whether the feature arises in circumstellar environments or in the surrounding molecular 
cloud, and constrain the volatility of the absorber. 



We would like to thank Dolores Walther for assistance with these observations, which 
were obtained during UKIRT Service Observing. We also appreciate useful conversations 
with Lou Allamandola, Christine Joblin, Scott Sandford, and Alan Tokunaga. 
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Table 1: Possible 3.25 Micron Feature Identifications 

Measured Measured Predicted A 
Source A (jum) T (K) at 80 K (/im) Ref. 



Mon R2/IRS-3 


3.249 ± 0.004 


80 


3.249 ± 0.004 


1 


Mon R2/IRS-3 


3.256 ± 0.003 


80 


3.256 ± 0.003 


2 


matrix-isolated coronene 


3.257 


4 


3.257 


3 


gas-phase coronene 


3.276 


698 


3.258 


4 


interstellar emission feature 


3.289 


1000 


3.260 


5 


hydrogenated amorphous carbon 


3.271 


300 


3.266 


6 


condensed coronene 


3.290 


788 


3.268 


4 


matrix-isolated pyrene 


3.268 


4 


3.269 


3 


quenched carbonaceous composite 


3.289 


300 


3.285 


7 



References — (1) Sellgren et al. (1994); (2) this paper; (3) Joblin et al. (1994); (4) Flickinger 
et al. (1991); (5) Tokunaga et al. (1991); (6) Biener et al. (1994); (7) Sakata et al. (1990). 

Note: The wavelength of these aromatic substances at a temperature of 80 K, appropriate 
for Mon R2/IRS-3 (Smith et al. 1989), was predicted from the measured wavelength and 
the temperature at which the wavelength was measured, using the temperature-dependent 
wavelength shifts measured by Joblin et al. (1995) for pyrene (for pyrene) or coronene (for 
all other substances). For the interstellar aromatic emission feature, we assumed a particle 
temperature of ~1000K (Sellgren, Werner, & Dinerstein 1983). 
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Figure Captions 



Figure 1 — New observations of the protostar Mon R2/IRS-3. Gaps in the data near 
3.30 /im and 3.32 /im are due to Pfund 5 in the standard star and strong telluric methane 
absorption, respectively. Top: the 3.16-3.55 /im spectrum (histogram) with a resolution of 
0.0033 fim (A/AA = 1000 at 3.25 /im). The units are flux density (Fx) in W cm -2 /tm" 1 vs. 
wavelength in microns. A third-order polynomial (solid curve) was fit to the observations, 
excluding 3.2-3.3 /im and 3.4-3.6 /im from the fit, to determine the continuum. Middle: 
the 3.16-3.55 /im optical depth (histogram) , compared to the sum of two Gaussians (solid 
curve), centered at 3.256 /im and 3.484 /im. The central wavelengths, widths, and optical 
depths of these two Gaussians were varied to produce the best fit to the data. Bottom: the 
3.16-3.55 /im optical depth (histogram) , compared to the profile of the aromatic interstellar 
emission feature (solid curve) in IRAS 21282+5050 (Nagata et al. 1988), after continuum 
subtraction (Tokunaga et al. 1991). The emission feature profile was first shifted to bluer 
wavelengths by 0.0294 /im to correct for temperature (see text and Table 1), and then 
scaled by the ratio of the average 3.17-3.28 /tm optical depth of Mon R2/IRS-3 to the 
average 3.17-3.28 /im feature profile of IRAS 21282+5050. 
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Figure 1 




Mon R2/IRS-3 and 
best-fit Gaussians 




Mon R2/IRS-3 and 
ISM emission feature 
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